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Abstract 

We discuss a new variant of the Eq inspired supersymmetric standard model (l^SSM) 
in which the two inert singlinos are exactly massless and the dark matter candidate 
has a dominant bino component. A successful relic density is achieved via a novel 
mechanism in which the bino scatters inelastically into heavier inert Higgsinos dur- 
ing the time of thermal freeze-out. The two massless inert singlinos contribute to 
the effective number of neutrino species at the time of Big Bang Nucleosynthesis, 
where the precise contribution depends on the mass of the Z' which keeps them 
in equilibrium. For example for mz> > 1300 GeV we find N e g ~ 3.2, where the 
smallness of the additional contribution is due to entropy dilution. We study a 
few benchmark points in the constrained -B^SSM with massless inert singlinos to 
illustrate this new scenario. 



1 Introduction 



The evidence for dark matter is now very strong. Initially proposed to allow an ex- 
planation of observed galactic rotation curves, we now also see its effects in the cosmic 
microwave background (CMB). The CMB baryon acoustic oscillation measurements from 
WMAP allow us to estimate the relative density of cold dark matter in the universe to 
be ficDM = (0.1099 ± 0.0062)//i 2 [1], where h is the reduced Hubble parameter h « 0.73. 

In supersymmetric (SUSY) models such as the minimal supersymmetric standard 
model (MSSM) [2] one typically imposes a Z 2 matter parity on the superpotential in 
order to remove the B — L violating terms allowed by the SM gauge symmetries. This 
is equivalent to so called i?-parity [3] , a Z 2 symmetry of the Lagrangian under which the 
charge of a physical state is related to its spin. The scalar and fermionic components 
of a chiral superfield have opposite i?-parity and the lightest supersymmetric (imparity 
odd) particle (LSP) is stable. In the such models one therefore predicts a new stable 
particle that could be a dark matter candidate, motivated for reasons not a priori related 
to dark matter. In general the LSP could be either a gravitino or a neutralino, depending 
for example of the nature of SUSY breaking which determines the typical value of the 
gravitino mass. 

In such theories a sub-weak-strength interacting neutralino is generally considered a 
good candidate for LSP dark matter [HE]. Since such a particle is self-charge-conjugate, 
its relic abundance is determined by thermal freeze-out, not by matter- antimatter asym- 
metry as in the case of baryons. At some point in the early universe the expansion 
rate of the universe would have become larger than the LSP's self-annihilation rate (or 
co-annihilation rate with other supersymmetric particles). At this point the number den- 
sity of LSPs would no longer be able to track its equilibrium value. It would remain 
much larger than the equilibrium value, only diluting due to Hubble expansion. An 
LSP dark matter candidate with a larger cross-section would be able to stay in equi- 
librium longer, meaning that there would be less of it in the universe today. Thermal 
neutralino dark matter has been widely studied in the MSSM [HI El El IS] and constrained 
(c)MSSM PHI EH EH EH EH- Successful dark matter may be realised if the LSP is the 
lightest neutralino, with various regions of parameter space corresponding to different 
dominant annihilation mechanisms. For example the bulk region corresponds to annilila- 
tion via t-channel slepton exchange, the focus region via t-channel chargino exchange and 
the funnel region via s-channel Higgs exchange. There are also other regions corresponding 
to co-annihilation with staus or stops. 

However the MSSM has certain shortcomings and it is possible that TeV scale SUSY 
is realised via a richer structure. It is therefore worthwhile considering alternative SUSY 
theories in which dark matter may be realised differently from in the MSSM, with the 
neutralino dark matter having a different composition and/or annihilating via different 
mechanisms. For example, the -E^SSM [TSJ EH EZ] (or the Exceptional Supersymmetric 
extension to the Standard Model (SM)) is a string theory inspired supersymmetric model 
based on an E s grand unification (GUT) group. The low energy gauge group contains 
an extra U(l), called U(1)n, under which the right-handed neutrinos that arise in the 
model are not charged. This means that the right-handed neutrinos may acquire large 
intermediate-scale Majorana masses, leading to a type-I see-saw mechanism to explain 
the small observed neutrino masses. This choice, that the low energy gauge group is 
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SM x U(1)n, defines the model. The Z7(1)jv g au g e symmetry is spontaneously broken at 
low energy by a SM-singlet field which we refer to as 5*3. This field radiatively acquires 
a vacuum expectation value (VEV) s naturally of order the soft SUSY breaking scale, 
leading to a Z'-boson with an induced mass of order the TeV scale. Automatic anomaly 
cancellation is ensured by allowing three complete 27 representations of E$ to survive down 
to the low energy scale. These three 27s contain the three generations of known matter, 
however they also contain the VEV acquiring Higgs doublets and SM-singlet. This means 
that there are two extra copies of the Higgs doublets and SM-singlet in the low energy 
particle spectrum. In the i? 6 SSM only one generation of Higgs doublets and singlets, 
defined to be the third, acquires the required VEVs and are called "active", namely H d3 , 
H u3 and S3. The other two generations, the first and second, of Higgs doublets H^a, H ua 
and SM-singlets S a , where a e 1,2, do not acquire VEVs and are called "inert". The 
inert generations have suppressed Yukawa couplings to SM matter, suppressing flavour 
changing neutral currents (FCNCs) and in turn explaining why they do not radiatively 
acquire VEVs. 

In the MSSM the neutralino mass matrix has the familiar 4x4 structure corresponding 
to the bino B, the neutral wino W 3 and two active^ neutral Higgsinos H® 3 and H® 3 . In the 
-E^SSM the neutralino mass matrix is greatly enlarged to a 12 x 12 structure including 
an additional four neutral inert Higgsinos H% a , H® a , one active singlino S 3 , two inert 
singlinos S a and an extra bino B' which is the SUSY partner of the Z' . It has been 
observed that the six inert Higgsinos and singlinos H® a , H® a and S a tend to decouple 
from the rest of the neutralino spectrum and it makes sense to consider their 6x6 mass 
matrix separately [T8]ll9]. Moreover the lighest inert neutralino state is predominetly inert 
singlino in nature and only acquires mass only via mass mixing with the inert Higgsinos 
proportional to the electroweak VEV v. It has a suppressed mass of order v 2 /s, where s 
is the SM-singlet VEV [18J. It is therefore natural to suppose that the LSP arises from 
the inert neutralino sector and is a superposition of the inert singlino and inert Higgsino 
components. It has been shown that it is possible to reproduce the observed relic density 
in this model by allowing annihilations via an s-channel Z-boson, while ensuring that the 
LSP, which must be similar in mass to the Z, having a mass of about 35-50 GeV, is not 
ruled out by LEP constraints on Z decays [TS]. It has recently been observed that, in 
this scenario, since the LSP always couples rather strongly to the SM-like Higgs boson, 
it leads to "invisible" Higgs decays at the LHC and for the same reason this dark matter 
scenario may be discovered or ruled out by dark matter direct detection experiments in 
the near future [lijffi 

In this paper we introduce a new scenario for dark matter in the E'gSSM in which the 
dark matter candidate is just the usual bino. At first sight having a bino dark matter 
candidate seems impossible since, as already discussed, the inert singlinos S a naturally 
have suppressed masses and it is very difficult to make them even as heavy as half the 
Z mass. It is clear that the singlinos will always be lighter than the bino, at least in 
the case of gaugino unification, as assumed in this paper. To overcome this we propose 

1 We shall refer to Higgsinos as being "active" or "inert" according to whether their scalar SUSY 
partners do or do not develop VEVs. 

2 Since this paper was submitted, the latest XENON100 direct detection results have been released [20] . 
These results severely challenge the previous dark matter scenario in jTSJ QI5] but are consistent with the 
bino dark matter scenario as discussed in this paper. 
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that the singlinos are exactly massless and decoupled from the bino, which is achieved in 
practice by setting their Yukawa couplings to zero. This is easy to do by introducing a 
discrete symmetry Zf under which the inert singlinos S a are odd and all other states are 
even, a scenario we refer to as the E^L\ SSM or EZSSM. In the EZSSM the inert singlinos 
S a will be denoted as a in order to emphasise their different (massless and decoupled) 
nature. The stable dark matter candidate (DMC) is then generally mostly bino and the 
observed dark matter relic density can be achieved via a novel scenario in which the bino 
inelastically scatters off of SM matter into heavier inert Higgsinos during the time of 
thermal freeze-out, keeping the bino in equilibrium long enough to give the desired relic 
abundance. Providing the inert Higgsinos are close in mass to the bino this is always 
possible to arrange, the only constraint being that the inert Higgsinos satisfy the LEP2 
constraint of being heavier than 100 GeV. This in turn implies that the bino must also 
be heavier than or close to 100 GeV. These constraints are easy to satisfy and, unlike the 
inert neutralino LSP dark matter scenario, we find that successful relic abundance can be 
achieved within a GUT-scale-constrained version of the model (the cEZSSM ), assuming 
a unified soft scalar mass tuq, soft gaugino mass M1/2 and soft trilinear mass Aq at the 
GUT scale [2D E21 E3] . 

It is worth noting that studies of the constrained (c)_E 6 SSM [211 E2J E3] have hitherto 
neglected to study the full 12 x 12 neutralino mass matrix, only considering the 6x6 
mass matrix of the MSSM augmented by the active singlino S3 and the extra bino B', as 
in the so called USSM |24J. Although the question of dark matter was addressed in the 
USSM, the requirement of successful relic abundance was not imposed on the ci^SSM 
[2Tj |22| 123] even though the latter analysis considered the same 6x6 mass matrix as in the 
USSM. Here we shall consider the c^SSM with the full 12 x 12 neutralino mass matrix, 
including both the USSM and inert neutralinos, under the assumption that the fermionic 
components of the inert SM-singlet superfields, the two inert singlinos, are forbidden to 
acquire mass by an extra Z 2 symmetry of the superpotential. In practice this reduces to 
a 10 x 10 neutralino mass matrix once the two massless inert singlinos are decoupled. 

In summary, the main result of this paper is that bino dark matter, with nearby 
inert Higgsinos and massless inert singlinos, provides a simple and successful picture of 
dark matter in the i^SSM consistent with GUT constrained soft parameters. We shall 
also consider the effect of the presence of the two massless inert singlinos in the EZSSM 
on the effective number of neutrinos contributing to the expansion rate of the universe 
prior to BBN, affecting 4 He production. Current fits to WMAP data [25] favour values 
greater than three, so the presence of additional contributions to the effective number of 
neutrinos is another interesting aspect of the EZSSM which we shall study. In practice 
we find that the additional number of effective neutrino species is less than two, due 
to entropy dilution, depending on the mass of the Z' which keeps the inert singlinos in 
equilibrium. 

The EZSSM is introduced in Section [2J Section [3] explores the neutralino and chargino 
sectors of the EZSSM whereas the details of the dark matter calculation are presented 
in Section |H N e g is defined and calculations of its value in the EZSSM are presented in 
Section |5j Some benchmark points are presented in Section [6] and the conclusions are 
summarised in Section [7] 
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2 The EZSSM 

The E e GUT group can be broken as follows: 

E 6 -> 50(10) x U{1)^ (1) 
->■ SU(h) x U{l) x x C/(l)v, (2) 
->■ SU(3) x ST/(2) x U(1) Y x C/(1) x x £7(1)^. (3) 

In the £ 6 SSM £ 6 is broken at the GUT scale directly to £77(3) x £77(2) x U(l) Y x C/"(1)jv, 
where 

17(1)* = cos(tf)tf(l) x + sm(tf)[/(l)^ (4) 

and tan($) = vTo" such that the right-handed neutrinos are chargeless. Three complete 
27 representations of Eq then survive down to low energy in order to ensure anomaly 
cancellation. These 27« (i £ {1,2,3}) decompose under the SU(5) x U(1)n subgroup as 
follows: 

27; ->■ (lO, l/v 7 !)) + (5,2/ V^j 

+ (d, . + (5, -2/Vio) + (l, 5/v 7 !)) + (l, 0) (5) 

The first two terms contain normal matter, whereas the final term, which is a singlet under 
the entire low energy gauge group, contains the right-handed neutrinos. The second-to- 
last term, which is charged only under U(1)n, contains the SM-singlet fields The third 
generation SM-singlet acquires a VEV (S 3 ) = s/V2 which, as we shall see, generates the 
effective /i-term and spontaneously breaks U(1)n leading to a mass for the Z'-boson. The 
remaining terms contain three generations of down- and up-type Higgs doublets H di and 
H ui as well as exotic coloured states D; and Di. Again only the third generation of Higgs 
doublets acquire VEVs {H%) = v d /V2 = v cos(/3) j\[2 and (77° 3 ) = v u /V2 = v sin(/3) /^2. 
The low energy gauge invariant superpotential can be written as follows: 

W = >v + >v 1 + h; 2 , (6) 

where 

Wo = \ijkSiHdjH u k + KijkSiDjDt + hfj k N^H U jLk 

+ hf jk u^H uj Q k + hf jk d c i H dj Q k + hf jk e1H dj L k , (7) 

Wi = g^D^Q. + g^D^ut, (8) 

m = gS^Djdt + gf^Djut + g^Q^D,. (9) 

It is now clear that the effective /i-parameter is given by /x = A 333 s/v^2 generating the 
term fj,H d3 H u3 in the superpotential. 

In order to suppress non-diagonal flavour transitions arising from the Higgs sector the 
superpotential obeys an approximate Z2 symmetry called Z^. Under this symmetry all 
superfields other than S3, H d3 and H u3 are odd. It is this approximate symmetry that 
distinguishes between the third generation and the inert generations of Higgs doublets and 
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SM-singlets, with the inert generations having suppressed couplings to matter and not 
radiatively acquiring VEVs. This approximate symmetry suppresses A^ couplings of the 
forms A a33 , A 3a3 , A 33a and A a/ g 7 , where a, /3, 7 e {1,2} i.e. labelling the inert generations 
only. Such an approximate Z^ symmetry, with a stable hierarchy of couplings, can be 
realised in .E 6 SSM flavour theories such as the one proposed by Howl et al. [2H]. The 
symmetry cannot be exact or else the lightest of the exotic coloured states would be 
absolutely stable which would contradict observation. 

Given this, an exact Z 2 symmetry must be imposed on the superpotential in order to 
avoid rapid proton decay. There are two ways to impose an appropriate Z 2 symmetry on 
W that leads to baryon and lepton number conservation. The first option is to impose 
a symmetry called Z 2 under which only the lepton superfields are odd. In this case the 
superpotential is equal to W0 + W1 and the model is called the ^SSM-I. U(1)b and U(1)l 
are symmetries of the superpotential if the exotic coloured states D and D are interpreted 
as diquarks and antidiquarks, with (B = ±2, L — 0). The second option is to impose 
a symmetry called Z^ under which both the lepton superfields and the exotic D and D 
superfields are odd. In this case the superpotential is equal to Wo + W 2 and the model is 
called the -B^SSM-II. U(1)b and U(1)l are symmetries of the superpotential if the exotic 
coloured states are interpreted as leptoquarks, with (B = =)=1, L = =)=1). 

It needs to be noted that the superpotential of the _E 6 SSM is already automatically 
invariant under the usual matter parity of the MSSM, provided that the exotic D and 
D superfields as well as the Higgs SM-singlet superfields are interpreted as being even 
under matter parity, along with the Higgs doublets. The B — L violating terms of the 
MSSM superpotential that matter parity is invoked to forbid are never present in the 
i^SSM superpotential of Eq. since they would violate the extra surviving U(l) n gauge 
symmetry contained in E e . We shall refer to this usual matter parity of the MSSM which 
automatically arises in this model as Z^ 7 . In the usual way it can be recast as -R-parity 
in terms of the superfield components, with the scalar components of the superfield being 
assigned the same parity as the superfield, and the fermionic components being assigned 
the opposite parity. As in the MSSM, the states which are odd under .R-parity are called 
the superpartners, with the lightest superpartner being absolutely stable. 

In the EZSSM the superpotential is also invariant under an additional exact Z 2 sym- 
metry called Zf. Under this symmetry only the two inert SM-singlet superfields S a 
are odd. The couplings of the forms \ a ij and K a ij are forbidden. This means that the 
fermionic components of S a , the inert singlinos a, are forbidden to have mass and do 
not mix with the other neutralinos. They only interact via their gauge couplings to the 
Z'-boson, which exist since they are charged under the extra U(1)n gauge symmetry. The 
extra Zf symmetry of the superpotential does not change the forms of the mass matrices 
of the the non-inert sectors of the model. In particular this means that all of the squark, 
slepton, gluino and non-inert Higgs scalar masses and mass matrices are the same as the 
ones given in [22] . The issue of Z-Z' mixing is also the same as discussed there. 

All of the exact and approximate Z 2 symmetries of the superpotential mentioned in 
this section are summarised in Table [U 

It is known that the model as presented thus far leads to gauge coupling unification at 
too high a scale, with the GUT scale typically being higher than the Planck scale. This 
issue can be solved by having the E e GUT group be broken to an intermediate group 
before being broken finally to SM x U(1)n as shown in [27]. In this paper, however, for 
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Table 1: The transformations of the superfields under the various Z2 symmetries of the superpotcntial 
that are mentioned in this paper. Z^ is an approximate flavour symmetry. Either Z2 or Zf is imposed 
in order to avoid rapid proton decay. Zf matter parity is already a symmetry of the i^SSM due to 
gauge symmetry. Zf is the extra symmetry which is imposed in the EZSSM, forcing the inert singlinos 
to be massless. 



simplicity, we implement the usual solution where the superpotential contains a bilinear 
term involving extra fields from incomplete 27' and 27' representations W = fi'H'H'. To 
some extent this solution reintroduces the \i problem, but // is not required to be related 
to the electroweak symmetry breaking (EWSB) scale and in order to observe satisfactory 
gauge coupling unification it is only required that jjt! < 100 TeV. 



3 The Neutralinos and Charginos of the EZSSM 

In the EZSSM the chargino mass matrix in the interaction basis 
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(10) 



(11) 



(12) 



(13) 



The top-left block is the MSSM chargino mass matrix whereas the bottom-right block 
contains mass terms for the extra inert Higgsino states. The Yukawa couplings in the 
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off-diagonal blocks are suppressed under the approximate and are therefore expected 
to be small, implying that there should not be much mixing between the MSSM and inert 
states. 

We define the term "neutralino" not to include the massless inert singlinos which have 
no Yukawa couplings involving them and are decoupled. The neutralino mass matrix M N 
in the interaction basis 
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(15) 



where once again a, (3 G {1,2}, indexing the inert generations. Qd = — -^=, Q u 
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and Qs = ^7= are the U(l) n charges of down-type Higgs doublets, up-type Higgs doublets 
and SM-singlets respectively and g[ is the GUT normalised U(1)n gauge coupling. M 1; 
M 2 and M[ are soft gaugino masses. Typically g[ ~ g\ all the way down to the low 
energy scale. If the soft gaugino masses are unified at the GUT scale (universal M 1 / 2 ) 
then we also have M[ pa M\ pa M 2 /2. The matrix as written neglects the small kinetic 
term mixing between B and B'. The elements left empty in this matrix and similar ones 
are implicitly taken to be zero. 

The Yukawa couplings in the off-diagonal blocks are suppressed under the approximate 
1j 2 . Given the smallness of these couplings, the inert neutralinos in the bottom- right block 
are pseudo-Dirac states with an approximately decoupled mass matrix 



s 
71 
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A322 
\ A321 



A322 
A312 



A312 
A311 



A321 \ 

^311 

/ 



in the basis ( H° d2 H° dl H° u2 H° ul ) 



They are approximately degenerate with the two inert chargino Dirac states. 

The top-left block contains the states of the MSSM supplemented by the third gener- 
ation singlino and the bino'. This is known as the USSM sector [28J. In the case where 
Mi pa M[ is small the lightest neutralino mass state will be mostly bino. The bino' will 
mix with the third generation singlino giving two mixed states with masses around Qsg'\S. 
As Mi pa M[ increases the bino mass will increase relative to both the third generation 
Higgsino mass fi and the inert Higgsino masses given approximately by the bi-unitary 
diagonalisation of — -^\z a ps. At the same time the state mostly containing the third 
generation singlino will have a decreasing mass as M[ increases relative to Qs9\S. 
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4 Dark Matter in the cEZSSM 



As discussed previously, due to the automatic matter parity of the model, there is a 
conserved i?-parity under which the charginos, neutralinos, inert singlinos a and exotic 
D and D fermions, along with the squarks and sleptons, are all i?-parity odd i.e. all 
of the fermions other than the quarks and leptons are .R-parity odd. We shall assume 
that the lightest neutralino N\ is the lightest of all of the -R-parity odd states, excluding 
the massless inert singlinos a. However Ni cannot decay into a via neutralino mixing 
since the inert singlinos are decoupled from the neutralino mass matrix. Furthermore the 
possible decay N% — > aa, allowed by the a-a-B' supersymmetric U{1)n gauge coupling, 
is forbidden if N\ is lighter than the inert SM-singlet scalars a. In fact, in this case, 
no kinematically viable final states exist that have the same quantum numbers as N±. 
Therefore N% is absolutely stable and is the DMC of the model. 

In the successful dark matter scenario presented in this section Ni has a dominant bino 
B component with at least one of the two pairs pseudo-Dirac inert Higgsinos expected 
to be close in mass, but somewhat heavier, in order to achieve the correct relic density. 
This is due to a novel scenario in which the DMC, approximately the bino, inelastically 
scatters off of SM matter into heavier inert Higgsinos during the time of thermal freeze- 
out, keeping it in equilibrium long enough to give a successful relic density. In this section 
we discuss in detail how this novel scenario comes about in this model. 



4.1 The Dark Matter Calculation 

Usually in supersymmetric models the evolution of the cosmological number density rii of 
a supersymmetric (.R-parity odd) particle i in the early universe can be expressed as 



j 

[( a 'xij v ix) ( n i n x ~ < q ^x) - WxjiVjx) (njnx - nfng) . (16) 

jfr x 

The first term accounts for Hubble expansion and the second term accounts for annihi- 
lations with other supersymmetric particles, including self-annihilations. The third term 
represents the decays of supersymmetric particles i into other supersymmetric species j 
as well decays of other supersymmetric species into species i. The final term represents 
the inelastic scattering of supersymmetric particles i off of R-parity even particles X into 
other supersymmetric species j and vice versa [29| 130]. 

Summing up these equations yields the somewhat simpler 

h = ^2hi = -3Hn - ^(vijVij) ( n i n j ~ n T n f) ■ ( 17 ) 

i i j 

It should be noted that, assuming all supersymmetric particles can decay into the DMC 
in a reasonable amount of time, after thermal freeze-out the relic number density of the 
DMC will subsequently becomes equal to n. 
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In our model the DMC is not the lightest .R-parity odd state (an inert singlino), but 
the lightest neutralino N\. We would like to use Eq. (fT6]) to describe the evolution of 
i?-parity odd states other than the inert singlinos, generically x- 111 this case we should 
also include in Eq. ( ITE)) processes involving a and b particles that change the number of 
X particles by one. Since such processes necessarily involve inert SM-singlet scalars a, it 
is valid to neglect these processes in the case where these inert SM-singlets have frozen 
out long before the freeze-out of dark matter. We will call this condition 1 and it should 
be satisfied given our assumption that the inert SM-singlet scalars are heavier than the 
DMC, since they only interact via the heavy Z'-boson. As we shall see, the value of n 
after the thermal freeze-out of Ni depends on annihilation cross-sections involving Ni and 
other .R-parity odd states close by in mass. As long as condition 1 is satisfied, meaning 
that we can neglect such annihilations that also have inert singlinos in the final state 
during thermal freeze-out, we can neglect the inert singlinos and use Eq. ( I17p to calculate 
the number density of i?-parity states other than inert singlinos. n will eventually be 
equal to the number density of DMCs after other x particles have decayed to N x . 

During thermal freeze-out the annihilation rates of the x particles become small com- 
pared to the expansion rate of the universe and their number densities become larger 
than their (non-relativistic) equilibrium values. The universe expands too fast for the 
number densities to track their equilibrium values. Let us assume however that these 
states inelastically scatter off of SM states X frequently enough that the ratios of the 
number densities of the x particles do maintain their equilibrium values during the time 
of thermal freeze-out. We shall call this condition 2 and assuming that it is satisfied we 
have 



n i = n l _ n i = n^_ 

eq : 

rii nf- n n eq 



which allows us to rewrite Eq. ( ITT)) as 

n = -3Hn - (av) (n 2 - n 2 cq ) , (19) 

where n eq = and 

eq eq 
Tl Tl 

M = ££<<w>-br-- (20) 



^eq 



Here we see that since heavier x particles would have smaller non-relativistic equilibrium 
number densities there would be fewer of them around during the dark matter's thermal 
freeze-out and annihilation cross-sections involving them would be less important. 



4.2 The cEZSSM 

In order to carry out the dark matter analysis in the constrained version of the model we 
have extended the RGE code used by Athron et al. [22] to include the Yukawa parameters 
and soft masses of the inert sector of the EZSSM. The inputs are k^j and A333 at the GUT 
scale, Xsap at the EWSB scale, s and tan(/3), as well as the known low energy Yukawa 
couplings and gauge couplings. Given these inputs and the RGEs the algorithm attempts 
to find points with GUT scale unified soft masses m , Mi/ 2 and A . The low energy 
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U(1)n gauge coupling g[ is set by requiring it to be equal to the other gauge couplings at 
the GUT scale. 

For consistent points in the -E^SSM the lightest non-inert (USSM sector) supersym- 
metric particle is typically bino dominated. For the cEZSSM we find the same thing. The 
masses of the inert Higgsino states depend on s and on the the Yukawa couplings \z a p 
and the in the cEZSSM the lightest neutralino can be either the bino dominated state or 
a pseudo-Dirac inert Higgsino dominated state. In the latter case we find that the DMC 
pseudo-Dirac inert Higgsino states co-annihilate with full-weak-strength interactions and 
lead to a too small dark matter relic density. In the former case the bino DMC normally 
annihilates too weakly and yields a too large dark matter relic density. If, however, there 
are inert Higgsino states close by in mass, they contribute significantly to (av) allowing 
for the observed amount of dark matter. This relies on condition 2 being satisfied i.e. the 
binos being up-scattered into inert Higgsinos with a large enough rate. 

Such points with consistent dark matter relic density can be found and three are 
presented in Section O Condition 1 is satisfied since the inert SM-singlet scalars are so 
much heavier than the DMC and the Z'-boson mass is so large compared to the regular 
Z-boson mass. Annihilation and scattering processes involving inert SM-singlets and 
singlinos must contain a virtual Z'-boson. 

To test condition 2 let us compare the rate for binos up-scattering into inert Higgsinos 
with the inert Higgsino co-annihilation rate. We shall label the mostly bino state N% and 
the lightest pseudo-Dirac inert Higgsino states N 2 and iV 3 . The dominant up-scattering 
diagrams are of the following form: 




We define Rzij couplings such that the Z-Ni-Nj coupling is equal to Rzij times the Z-v-v 
coupling. We have 

R Zij = £ N ? N ?~ £ W» ( 21 ) 

D=3,7,9 £7=4,8,10 

where N°" is the neutralino mixing matrix element corresponding to mass eigenstate i and 
interaction state a. D and U index the down- and up-type Higgsino interaction states 
respectively. For the pseudo-Dirac inert Higgsino states we have 

m 3 « -m 2 and R Z 23 ~ 1, (22) 
allowing for full-weak-strength co-annihilations of the following form: 
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The ratio of the rate for the mostly bino state up-scattering into the mostly inert 
Higgsino state to the inert Higgsino co-annihilation rate is given approximately by 



T 



(cT23^23)^2 q % 




(23) 



To give an idea of the size of this ratio, if the SM particle X is relativistic and mi ~ m 2 « m 3 



where x = \m\\/T and T is the temperature. This ratio is expected to be large because 
of the overwhelming abundance of the relativistic SM particle X, but it also depends on 
Rzn- The value of Rz\2 depends on the Z,f -breaking couplings that mix the top- left 
block of the neutralino mass matrix in Eq. ( I15p . the USSM states including the bino, with 
the inert Higgsino states in the bottom-right block. Since this symmetry is not exact we 
expect these couplings to be large enough such that we can still assume T ^> 1. Explicit 
examples of this parameter are included in Table [3] in Section [61 

With the two conditions satisfied we use micrOMEGAs [31 J to calculate the dark matter 
relic density for low energy spectra consistent with the GUT-scale-constrained scenario. 
The CalcHEP model files for the EZSSM are produced using LanHEP [32]. The observed 
relic density of dark matter can arise in this model and examples are shown in Section [61 
The most critical factor is the mass splitting between the bino and the lightest inert 
Higgsinos. Too large and there would not be enough inert Higgsinos remaining at the 
time of the bino's thermal freeze-out to have a significant enough effect. Too small and 
(crv) would be dominated by inert Higgsino co-annihilations leading to a too small dark 
matter relic density. 

Since in this scenario the DMC is predominantly bino, the spin-independent DMC- 
nucleon cross-section cr SI is not expected to be in the range that direct detection exper- 
iments are currently sensitive too. The spin-independent cross-section of a pure bino is 
suppressed by the squark masses, but is also sensitive to the squark mixing angles [35] . 
For each flavour the cross-section vanishes for zero squark mixing. Since in practice the 
DMC will also have non-zero (but small) active Higgsino components, there are also con- 
tributions to o"si from t-channel active Higgs scalar exchange, via the bino-Higgs-Higgsino 
supersymmetric gauge coupling. These contributions, though dominant, are quite small 



then 



T 





3/2 



(24) 
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due to the overwhelming bino nature of the DMC Estimates of crsi, using the same pro- 
ton fd, f u and f s parameters used by Gogoladze et al. [31], are included in Table [3] in 
Section 

5 The Inert Singlinos and Their Contribution to the 
Effective Number of Neutrinos prior to BBN 

In the standard theory of BBN, which happens long after the thermal freeze-out of 
dark matter, the resultant primordial abundances of the light elements depend on two 
parameters — the effective number of neutrinos contributing to the expansion rate of the 
radiation dominated universe N e g and the baryon-to-photon ratio r\. 

Whilst the primordial abundance of 4 He is not the most sensitive measure of 77, it 
is much more sensitive to N e g than the other light element abundances. This is because 
prior to nucleosynthesis when the equilibrium photon temperature is of order 0.1 MeV the 
number of neutrons remaining, virtually all of which are subsequently incorporated into 
4 He nuclei, is sensitive to the expansion rate of the universe, which depends on N e s. The 
greater the expansion rate, the less time there is for charged current weak interactions to 
convert neutrons into protons. 

The analysis by Izotov et al. [35] using the more recent neutron lifetime measurement 
by Serebrov et al. [36] gives N c s = 3.80^q;7q at 2-sigma, implying a more-than-2-sigma 
tension between the measured 4 He abundance and the Standard Model prediction for N e g 
(about 3). Although Aver et al. [37J suggest that these errors may be larger, similar results 
are also obtained for the effective number of neutrinos contributing to the expansion rate 
of the universe from fits to WMAP data [25J . 

In the EZSSM the two massless inert singlinos would have decoupled from equilibrium 
at an earlier time than the light neutrinos, but nevertheless would have contributed to 
N c ff. Exactly when the inert singlinos would have decoupled from equilibrium with the 
photon depends on the mass of the Z'-boson which determines the strength of an effective 
Fermi-like 4-point interaction vertex that would have been responsible for keeping the inert 
singlinos in equilibrium. The various values for N e g that can be achieved in this model 
all fit the data better than the SM value. 

The implications of extra neutrino-like particles present in the early universe have 
long been studied and the methods used in following analysis rely on relatively simple 
physics [3S]- The cosmological energy density of a relativistic boson, fermion i with number 
of degrees of freedom g l and temperature T l is given by 

p i = (i,w^on 4 . (25) 

The total radiation dominated energy density is then defined to be 

p = g^T\ (26) 

where T = T 7 is the photon temperature. The effective number of degrees of freedom g 
takes into account the factor of 7 /s for fermions and also takes into account the fact that 
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some species no longer in equilibrium with the photon may have a different temperature. 
In this radiation dominated universe the expansion rate is then given by 



H 2 



8ttG 



(27) 



3 

f 



P 



off 



4vr 3 



(28) 




45 



(29) 



where we define the constant k\ for future convenience. 

The effective number of degrees of freedom contributing to the expansion rate of the 
universe during the run-up to nucleosynthesis is defined to be 



Here g 1 = 2 is the number of degrees of freedom of the photon and g v = 2 is the 
number of degrees of freedom of a light neutrino. The three SM neutrinos are expected 
to decouple from equilibrium with the photon at a temperature above the electron mass 
whereas nucleosynthesis does not happen until the temperature is below the electron 
mass. When the photon/electron temperature is around the electron mass the electrons 
and positrons effectively disappear from the universal. Their disappearance heats the 
photons to a higher temperature then they would otherwise have had, but the neutrinos, 
having already decoupled, would continue to cool at the full rate dictated by Hubble 
expansion. Because of the neutrinos' lower temperature at nucleosynthesis they would 
contribute less to per degree of freedom. Eq. f[3"Uj) is defined such that in the SM 
iV e g = 3, for the three neutrinos decoupling above the electron mass (as we shall see). 
Extra particles, such as the EZSSM inert singlinos, decoupling above the muon mass 
would have even lower temperatures at the time of nucleosynthesis and would therefore 
contribute to g$ s even less than light neutrinos per degree of freedom. 



In the cEZSSM there is a typical scenario in which the massless inert singlinos a decouple 
at a temperature above the colour transition temperature (when the effective degrees of 
freedom are quarks and gluons rather than mesons) and above the strange quark mass, 
but below the charm quark mass. This has to do with the strength of the interactions that 
keep the inert singlinos in equilibrium which depend heavily on the mass of the Z'-boson 
mass. If the inert singlinos do decouple in this range, this leads to a definite prediction 
for N e g. We shall explain why the inert singlinos typically decouple in this temperature 
range in the next subsection. For now we derive the value of N e s in this scenario as an 
example. 

We shall use the subscript to denote quantities at some temperature To below the 
electron mass and the subscript e to denote quantities at some temperature T e above 
the electron mass and where all light neutrino species are still in equilibrium. We shall 

3 A much smaller number of electrons remains due to the small lcpton number asymmetry. 



.eff 



g-y + V8^iV eff (4/ n ) 4 /3. 
2 + V4iV eff (Vn) 4/3 



(30) 
(31) 



5.1 The Calculation of N e g 
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use the subscript s to denote quantities at some still higher temperature T s above the 
colour transition and the strange quark mass and where the inert singlinos are still in 
equilibrium. 

At T s the effective number of degrees of degrees of freedom contributing to the expan- 
sion rate is 

gf = gl + g 9 + 7/ 8 ( g e + g » + g u + g d + g s + 3g v + 2g * ) 

= 2 + 16 + 7/s(4 + 4 + 12 + 12 + 12 + 6 + 4) = (32) 



and at T e it becomes 



and at T it becomes 



gf = 2 + V8(6 + 4[^-) 1 (33) 



rpv \ 4 / rpc \ ^ 

-L n \ . I J- n 



9 f = 2 + V8 ^^j + 4^j j, (34) 

taking into account that the neutrinos and inert singlinos now have different temperatures. 

The entropy within a given volume V due to a relativistic boson, fermion % with number 
of degrees of freedom g % is given by 

2tt 2 

S* = (1, 7 A)^— (T) V. (35) 

Since we are assuming that the inert singlinos decouple before the strange quark threshold, 
in going from T s to T e we conserve the entropy in the co-moving volume separately for 
the inert singlinos and for everything else. Specifically for the inert singlinos 

T!V S = (Tf) 3 K (36) 

and for everything else 

\9 1 + 9 g + 7 M9 e + (f + 9 v + 9 d + g 8 + Z9 v )]T?V 9 = [^ + 7 8 (/ + 3<f)]T e V e 

613/4TJV, = 103/ 4 T e V e . (37) 

This allows us to write 

TX _ fT!\ 3 _ 103/4 _ 43 



T e 3 K \T e J 613/4 247 

In going from T e to T we conserve the entropy separately for the neutrinos, for the inert 
singlinos again, and for everything else 

W + y 8 g e ]TX = giT*V , (39) 
T e V e = (T») 3 V , (40) 
(Tf)V e = (7*)V . (41) 
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This gives us 



To J g^ + 7 /sg e 11 



(42) 



and 



HY 43 a 1 43 4 f43 x 

T J 247 g^ + 7 /sg e 24711' 1 J 

In this case the effective number of neutrinos contributing to the expansion rate prior to 
nucleosynthesis (at T ) is then 

/ 43 y/ 3 

NeS = 3 + 2 v247; ^ 3 ' 194 ' (44) 

5.2 The Inert Singlino Decoupling Temperature 

The light neutrinos are kept in equilibrium via their electroweak interactions. For all the 
light neutrinos there are the following tree-level diagrams: 




For the electron neutrinos there is also the following additional diagram: 




We express the cross-section for processes relevant for keeping muon and r neutrinos in 
equilibrium as 

T 2 (5M2 4 

where k 2 is another constant and 

X 4 = ((-y 2 + 4,)(y 2 )) 2 + 1/4^ » 0.031. (46) 
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Note that using the GUT normalised U(l)y gauge coupling g\ we have 

92_ = /5 gi 
cw V 3 sw 

The cross-section for electron neutrinos with their extra diagram is then 



(47) 



T 2 (5/s) 2 fl 4 

M = k 2 ^ { -^Y\ (48) 
m z s w 

where 

Y 4 = ((i/2 + s 2 w )(^)) 2 + 1/4^ « 0.147. (49) 

We express the number densities of all Weyl fermions still in equilibrium with the photon 

as 

n e L = n e R = n » L = n m = n v e = n v„ = tfr = ksT 3 ^ 

and the expansion rate is given by 

H = h^gfT 2 . (51) 
The neutrino decoupling temperature T v can then be approximated by 

(a v v)n v = H (52) 

=► (T^f = Ky^fm% J ^l i , (53) 

(T*) 3 = K^gfm^^-^ (54) 

with K = /ci//c2^3- A more detailed calculation finds that in the SM (with only neutrinos, 
electrons and photons contributing to gl ) T v >" Vt 3.7 MeV and T Ue m 2.4 MeV, the 
muon and r neutrinos decoupling earlier. 

At temperatures above the strange quark mass the processes relevant for keeping the 
inert singlinos in equilibrium are as follows: 



e L ,e R ,fj, L ,fj, R ,u e ,u^,u T , 
UL,u R ,d L ,d R ,s L ,s R 



ul, u r , d L , d R , s L , s R 




The part of the Z' current illustrating the relevant U(1) N charges is 

/ (2)L \ 

(1) e* 

r z , = (l eR q u R d R $1 

(2) d R 
V (5)* J 



(55) 
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and the total cross-section taking into account all of these diagrams is then (neglecting 
the small Z-Z' mixing) 

T 2 Z 4 

(W) = h^gf-^, (56) 

where 

Z 4 = (5) 2 [2(2) 2 + 2(1) 2 + 3(1) 2 + 3(1) 2 + 6(1) 2 + 6(2) 2 + 3(2) 2 ] = 1450, (57) 
leading to an approximate singlino decoupling temperature of 

1 (40)2 



(T a ) 3 = Ky/jfm^^j- (5* 



T*y j g f (m z ,\ A m 2 m 2 Y 4 (59) 



T"'J ygf\m z J Z*' 

The only unknown variable here affecting the inert singlino decoupling temperature is 
then the Z 1 mass m,z>- Rearranging we find 



\ 3/4 



^ W m ^6^MevJ ■ (60) 



5.3 N eS in the EZSSM 

We now check which values of mz> are consistent with our assumption that the inert 
singlinos decouple at a temperature between the strange and charm quark masses. For 
T u < m c we find that we require mz> < 4700 GeV. For mz> ~ 1000 GeV the situa- 
tion is slightly more complicated. Firstly the temperature of the QCD phase transition 
is not accurately known and secondly the effective number of degrees of freedom is de- 
creased by so much after the QCD phase transition that even if the inert singlinos were 
decoupled beforehand the universe may be expanding slowly enough afterwards that they 
could come back into equilibrium. After checking a range of scenarios we find that for 
1300 GeV < m z , < 4700 GeV our value of N cS = 3.194 is valid. For m z > < 950 GeV 
the inert singlinos decouple at a temperature above the muon mass, but below the pion 
mass leading to a larger prediction of N c g = 4.373. The current experimental limit in the 
EZSSM is mz' > 892 GeV [39], so at the time of writing it is possible that the Z'-boson is 
light enough to predict N e g = 4.373. For Z' masses in between these ranges the value of 
N cS depends on the details of the QCD phase transition, but is somewhere between these 
predictions. For inert singlinos decoupling above the pion mass, but after the QCD phase 
transition we have N e g = 4.065. All of these values are within the 2-sigma measured 
range iV e fj = 3.80l^7o an< ^ cl° ser to the central value than the SM result N e s = 3. 



6 Benchmark Points 

In the following tables we present three benchmark points in the cEZSSM. For all three 
points we fix A 322 = 0.1 and A 32 i = A 3 i 2 = 0.0001 at the EWSB scale. For the Z^- 
breaking couplings we also fix A 332 = A 323 = 0.0 1 2 and A 33i = A 3i3 = 0.005 at the EWSB 



17 



Benchmark 


1 


2 


3 




j / /-)\ 

tan(p) 


30 


10 


3 




s [TeV] 


5 


4.4 


5.5 


^333 


@ GUT scale 


-0.3 


-0.37 


-0.4 


A322 @ 


EWSB scale 


0.1 


0.1 


0.1 


A311 @ 


EWSB scale 


0.0293 


0.0403 


0.0399 


^3ii 


@ GUT scale 


0.18 


0.18 


0.23 




M 1/2 [GeV] 


590 


725 


908 




m [GeV] 


1533 


454 


1037 




A [GeV] 


1375 


1002 


413 



Table 2: The GUT scale parameters of the three benchmark points. 





Bene! 


imark 


-1 

1 


O 

2. 




3 






[GeV] 


-lUoD. ( 


-1 loy.o 


-14U0.0 


A 3 22s/v / 2 


[GeV] 


353.55 


331.13 


388.91 


A311 


s/V2 


[GeV] 


103.59 


125.38 


155.17 


N, 


mass 


[GeV] 


94.07 


114.49 


143.50 


N 2 


mass 


[GeV] 


-105.12 


-126.45 


-156.57 


iV 3 


mass 


[GeV] 


105.14 


126.47 


156.62 


V 4 


mass 


[GeV] 


167.05 


203.19 


255.47 


N 5 


mass 


[GeV] 


-353.77 


-311.29 


-389.12 


N 6 


mass 


[GeV] 


353.78 


311.30 


389.13 


N 7 


mass 


[GeV] 


-1092.5 


-1194.5 


1409.6 


N s 


mass 


[GeV] 


1093.3 


1194.8 


-1411.2 


Ng 


mass 


[GeV] 


-1803.2 


-1572.3 


-1964.7 


N w 


mass 


[GeV] 


1899.7 


1688.7 


2109.9 


Ci 


mass 


[GeV] 


105.04 


126.41 


156.52 


C 2 


mass 


[GeV] 


167.05 


203.19 


255.46 


£3 


mass 


[GeV] 


353.78 


311.30 


389.13 




mass 


[GeV] 


-1094.4 


-1196.1 


-1411.3 




m z > 


[GeV] 


1850.4 


1628.4 


2035.4 






N cS 


3.194 


3.194 


3.194 






0.112 


0.107 


0.102 






T 


1.1 x 10 8 


2.3 x 10 s 


2.3 x 10 8 






[cm 2 ] 


4.9 x 10~ 48 


2.5 x 10" 48 


1.2 x 10" 48 



Table 3: The low energy neutralino and chargino masses, and associated parameters. The dark matter 
candidate is the lightest neutralino Ni, which is predominantly bino. There is a nearby pair of inert 
neutral Higgsinos N2, N3 and a chargino C\ into which N\ inelastically scatters during freeze-out, resulting 
in the correct relic density Ocdm'* 2 shown. The predicted values of mz 1 and N e g are also shown, as is 
the spin- independent Ni direct detection cross-section 031 ■ 
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Benchmark 


1 

X 


9 




il\ Illclbb 


[KjC V J 


1 99 9 


114 6 


11J.O 




1 VJt V J 




Q87 1 


1 ^99 


/l 1 > 1 Q G G 


[vjL V J 


1890 


1664 


2080 


/- T — 1— T V"l , \ t 1 t 1 

xl Illclbb 


[LrC V J 


21 06 


1 3Q6 


1 67^ 

XVJ ( tj 


y l 111CI00 


GeV 


91 rn 


1 ^Q^ 




lHS 2 i ,n S\ 




1547 


518 


1084 


m H d2 , m H dl 


1 VJfC V J 


1 567 


61 1 


1156 

X X ' .7 V 


m Hu2 ,m Hul 


PpV 
[KjC V J 


1561 


599 

U U 


1146 




[>aC V J 


1 

1100 


ouo 


1 7Q/L 




TPpV 
[KjtC V J 


1 443 

X 11 ■ J 


493 


1 77^ 




rnpv 

[VjL V j 


9864 


9^91 


^06^ 




[uC V J 




zoio 


ouoz 




[GeV] 


1 1 00 
1 LA A 


roc; q 


1 1 1 n 




[GeV] 


1817 


1774 


1707 




[GeV] 


1470 


1069 


1546 




[GeV] 


1838 


1224 


1761 




[GeV] 


1434 


1009 


1512 


m Sl ,m di 


[GeV] 


1840 


1226 


1763 


m l 2 


[GeV] 


1748 


1265 


1818 




[GeV] 


1907 


1278 


1820 




[GeV] 


1500 


718.8 


1259 




[GeV] 


1655 


731.3 


1261 


m f2 


[GeV] 


1708 


949.2 


1473 




[GeV] 


1775 


952.8 


1474 




[GeV] 


1705 


945.6 


1472 




[GeV] 


1774 


949.5 


1472 


rrig 


[GeV] 


541.3 


626.9 


787.7 



Table 4: The remaining particle spectrum in the standard notation. 



19 



scale. At the GUT scale we fix K333 = K322 = K311 and = for i ^ j. The lightest 
(SM-like) Higgs mass is calculated to second loop order. 

We have chosen three points with quite different values of tan(/3) — 30, 10 and 3. This 
illustrates the fact that tan(/3) can be quite low in this model since the SM-like Higgs 
mass is not constrained to be less than mz\ cos(2/3)| at tree-level as it is in the MSSM. 

The mass of the bino DMC N\ is not directly constrained to be above above 100 GeV. 
However, the lightest pseudo-Dirac inert Higgsino neutralinos N2 and V3 are almost de- 
generate with the lightest inert Higgsino chargino C\ and therefore these are constrained 
to heavier than 100 GeV in order to be consistent with LEP constraints. Furthermore the 
thermal relic DM scenario outlined in Section H] requires N 2 and iV 3 not to be too much 
more massive than N\. In practice the N% is predominantly bino and its mass cannot be 
much less than 100 GeV. In Benchmark 1, for example, it is 94 GeV. 

Requiring such values for the low energy bino mass M\ and requiring consistent elec- 
troweak symmetry breaking in practice means that the SM-singlet VEV s cannot be too 
low. This in turn means that the Z' mass is always quite a bit above the experimental 
limits, more than about 1.5 TeV. In these benchmarks from the constrained scenario the 
effective number of neutrinos contributing to the expansion rate of the universe prior to 
BBN iVeff therefore takes on the lower value calculated in Section [5] — around 3.2. This is 
more consistent with data than the SM prediction. 

In all benchmark points N4 and C2 are predominantly wino. iV 5 , iV 6 and C3 are 
predominantly made up of the rest of the inert Higgsinos states, with masses around 
A322.s/v / 2, whereas N?, N% and C4 are predominantly made up of the active Higgsinos 
states, with masses around \i. Ng and N w are mostly superpositions of the active singlino 
and bino'. 

The fact that T ^> 1 indicates that the inert Higgsino components in the predomi- 
nantly bino state N±, though small, are large enough such that processes involving N\ 
up-scattering off of a SM particle into N2 happen overwhelmingly more often than neu- 
tralino (co-)annihilation processes. In this way the ratios of the number densities of these 
particles are able to maintain their equilibrium values. 

The spin-independent DMC-nucleon cross-section cr SI , as estimated using the results 
of Choi et al. [33], is quite small for these benchmarks, and is not currently detectable by 
direct detection experiments. This is due to the predominantly bino nature of the DMC, 
and the large squark masses. 

7 Conclusions 

The question of dark matter in the -B^SSM illustrates the interesting diversity of possi- 
bilities that that can arise once one goes beyond the MSSM. The difficulty in making the 
inert singlinos predicted by this model much heavier than 50 GeV makes them natural 
dark matter candidates, but also led to a very tightly constrained scenario in which the 
inert LSP (essentially a mixture of the inert singlinos and inert Higgsinos) is now severely 
challenged by the most recent XENON100 analysis of 100.9 days of data. Moreover the 
tightly constrained parameter space makes it practically impossible for such a scenario to 
be consistent with having universal soft mass parameters. 

In this paper we have discussed a new variant of the _E 6 SSM, called the EZSSM, 
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which involves a novel scenario for Dark Matter in which the dark matter candidate is 
predominantly bino with a mass close to or above 100 GeV which is fully consistent with 
XENON100. A successful relic density is achieved via its inelastic scattering into nearby 
heavier inert Higgsinos during the time of thermal freeze-out. The model also predicts 
two massless inert singlinos which contribute to the number of effective neutrino species 
at the time of Big Bang Nucleosynthesis, depending on the mass of the Z'-boson which 
keeps them in equilibrium. For example for mz> > 1300 GeV we find N e s ~ 3.2. 

We have studied a few benchmark points in the constrained EZSSM with massless 
inert singlinos to illustrate this new scenario. The benchmark points show that it is easy 
to find consistent points which satisfy the correct relic abundance as well as all other 
phenomenological constraints. The points also show that the typical Z' mass is expected 
to be around 2 TeV, with the gluino having a mass around 500-800 GeV and squarks and 
sleptons typically having masses around 1-2 TeV. The direct detection spin-independent 
cross-sections agi ~ few x 10~ 48 cm 2 are well below current sensitivities. 
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